A multiband envelope function model for quantum transport in a tunneling diode 



Omar Morandi^B and Michele Modugno 2 

1 Dipartimento di Elettronica e Telecomunicazioni, Universita di Firenze, 
Via S. Marta 3, 50139 Firenze, Italy 
2 LENS, Dipartimento di Fisica and INFM, Universita di Firenze, 
Via Nello Carrara 1, 50019 Sesto Fiorentino, Italy 
(Dated: February 6, 2008) 

We present a simple model for electron transport in semiconductor devices that exhibit tunneling 
between the conduction and valence bands. The model is derived within the usual Bloch-Wannier 
formalism by a fc-expansion, and is formulated in terms of a set of coupled equations for the electron 
envelope functions. Its connection with other models present in literature is discussed. As an 
application we consider the case of a Resonant Interband Tunneling Diode, demonstrating the 
ability of the model to reproduce the expected behaviour of the current as a function of the applied 
voltage. 

PACS numbers: 85.30.Mn, 73.40.-c 



I. INTRODUCTION 

In recent years there has been a growing interest for 
semiconductor devices characterized by tunneling effects 
between different bands, as the Resonant Interband Tun- 
neling Diode (RITD) p|. This kind of diode belongs to 
the class of heterostructures that show a negative differ- 
ential resistance in a certain range of the applied volt- 
age, like the widely employed RTD (Resonant Tunneling 
Diode) @, OJ. However, differently from the latter where 
the electronic current flows within a single band, the re- 
markable feature of a RITD is the possibility to achieve 
a sharp coupling between "conduction" and "valence" 
states, allowing an interband current which becomes the 
main transport phenomena in the resonant region. 

The description of electron transport in such quantum 
devices hence requires multiband models capable to ac- 
count for tunneling mechanisms between different bands 
induced by the heterostructure design and the applied 
external bias. 

In the literature, different methods are currently em- 
ployed for characterizing the band structures and the 
electronic or optical properties of these heterostructures, 
such as envelope functions methods based on the effective 
mass theory [j, d, 0] , tight-binding 0, [1] and pseudopo- 
tential @ methods. In addition, various mathematical 
tools are employed to exploit the multiband quantum dy- 
namics underlying the previous models: the Schrodinger- 
like models [101 ] . the nonequilibrium Green's function 
[ill, 0, the Wigner function approach [i| H 0, 
and recently the hydrodynamics multiband formalisms 

[ana. 

All of these methods rely on some common approxima- 
tions to account for the effects of a non-uniform band pro- 
file on the electron dynamics. In particular, in the usual 
"k • P" approach [18|] , one starts by defining the Hamil- 
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tonian matrix of the bulk (in fc-space), and then allows 
the physical parameters (typically the band eigenvalues 
or the Luttinger-Kohn parameters [l9|) to have some x- 
dependence in order to describe the position-dependent 
properties of the heterostructure. In this approach care 
must be taken to preserve the self-jointncss of Hamilto- 
nian matrix, so appropriate quantization rules are needed 
[20l | . In this way the electrical fields arising from the 
band edge offset among different layers are not included 
from the beginning in the derivation of the model, but 
appear only at the macroscopic level (i.e. at the level 
of envelope functions). Indeed, in the previous approxi- 
mation technique, the ^-dependence of the unperturbed 
Hamiltonian matrix elements, generates a "mean effec- 
tive electric field" acting on the envelope functions, which 
is not present at the microscopic level. 

A different approach has been proposed in [2l[ were a 
local "modified Wannier basis" is chosen to include the 
inhomogeneity directly into the basis elements. Unfortu- 
nately, in this case the equations of motion of the enve- 
lope functions depend on the change of the Bloch func- 
tions across the interfaces (that are implicitly neglected 
in the previous procedure) and such an evaluation can 
result in a very difficult task. 

In this paper we introduce a different strategy, describ- 
ing the band edge offsets by means of external potentials 
applied to the bulk structure. This allows us to treat 
on the same footing both the electrostatic potential gen- 
erated by the charge distribution in the device and the 
heterostructure design of band edges, in order to high- 
light the role played by the heterostructure potential in 
the interband tunneling process. 

Within this framework we derive a hierarchy of multi- 
band models obtained by means of a fc-expansion, were 
the momentum fc plays the role of asymptotic parameter 
as in the usual "k-P" approach. The starting point is the 
single electron Bloch representation, that here we con- 
sider for simplicity for the case of non degenerate bands 
and constant band gaps. Then, after the fc-expansion, 
the electron wavefunction is projected on the Wannier 
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basis, yielding a set of coupled Schrodinger equations for 
the electron envelope functions in definite energy bands. 

These equations share some similarities with those of 
the well known Kane [13] and Luttinger-Kohn (LK) [III] 
models. However, a key difference is the choice of the 
basis elements. Indeed, since in a uniform crystal the 
Wannier functions of a given energy band are related to 
the Bloch functions of the same band by a unitary trans- 
formation, this allows us to give a simple physical insight 
to the envelope functions. Differently, since the Kane 
model arises form a unitary transformation of only the 
periodic part of the Bloch functions, the generic element 
of the Kane basis is nondiagonal in the Bloch band in- 
dex n, and envelope functions related to different "band" 
indices turn out to be coupled even in absence of any 
applied potential, therefore lacking of a direct physical 
interpretation. 

The LK model instead is a multiband effective mass 
model obtained from the latter by an additional quasi- 
unitary transformation that removes the spurious inter- 
band coupling to first order in k. However, since the LK 
approach is devoted to describe intraband effects, the 
coupling due to the external field is generally neglected. 

As an application of the present approach we consider 
the case of a two-band RITD, showing that the model is 
able to reproduce the expected behaviour of the current 
as a function of the applied voltage. 

The paper is organized as follows. In the next section 
we discuss the derivation of the model and the approxi- 
mations employed; then in Sec JIIII we analyze differences 
and analogies with the Kane [22l| and Luttinger-Kohn [ijj 
models. Finally, in Sec. IIVI we work out explicitly the 
case of a RITD, investigating its current-voltage charac- 
teristic curve. 



II. DERIVATION OF THE MODEL 

Let us consider an electron of mass m immersed in a 
crystal lattice described by the periodic potential Vl, in 
the presence of an additional external potential U that 
will be treated as a perturbation. The evolution of the 
electron wavefunction ^(x, t) is given by the solution of 
the Schrodinger equation 



2m 



v 2 + y L (x) + c/(x) 



*(x,t). (1) 



The eigenfunction of the unperturbed Hamiltonian 
#o = —(h 2 /2m)V 2 + Vl are Bloch functions Vn(k,x) 
(see e.g. |18j ) 



H tp n (k,x) = £ n (k)V>„(k,x) 



(2) 



and form a complete set with the orthonormality condi- 
tion 



f C(k,x)Vv(k',x) = 5(k-k% 

J X 



n being the band index and k the electron quasimomen- 
tum. Eq. (flj can be transformed in momentum space 
by means of standard textbook methods [23|, [13] that we 
review here below for completeness. According to the 
Bloch theorem the Bloch functions VVi(k, x) can be writ- 
ten as 



V>n(k,x) = e zK ' x ii„(k,x) = (x|n,k) 



(4) 



where the functions it„(k, x) have the same periodicity 
of the lattice potential and are normalized according to 



(27T) 



■ / u*(k,x)u„/(k,x 

J cell 



) — $nn' ■ 



" J cell 

A generic solution of Eq. (fT]) can be expanded as 

<&(X, t)=J2j k Vn{K t)^„(k, X) 



(5) 



(6) 



where k runs over the first Brillouin zone; then the ex- 
pansion coefficients satisfy the following equation (here- 
inafter we omit the time dependence to simplify the no- 
tation) 

ihdtM*) = E n (k)<p n (k) + f {n,k\U\n',k')<p n ,(k'). 

n' 

(7) 

By exploiting the periodicity of the u n (k, x) functions 
the expectation value of the external potential U can be 
rewritten as 

(n,k\U\n',k') = f e i ( k '- k > x <(k,x)ii„,(k',x)[/(x) 

= ^A(n,<k,k') / e l ( k '- k - K ')' x [/(x) 
1 J* 

= (2n) 3 Y / B iU(k' -k-K t ) (8) 
1 

K.i being a reciprocal lattice vector, and U the Fourier 
transform of U. 

At this point, following [23| . we assume the potential U 
to be nearly constant over a single lattice cell, so that only 
the zero momentum Fourier component give a relevant 
contribution 

(n, k\U\n', k') ~ (2ir) 3 B U{k' - k) (9) 

where Bg can be expressed as 

B o = / <(k,x)u„/(k',x) = (w n kK' k') ( 10 ) 

" J cell 

f2 being the volume of a single cell. 

Let us now evaluate explicitly the coefficients Bq, by 
considering separately the case n — vl and n ^ n' . In 
the former case it is easy to show from Eq. that 



(3) 



Bo(n,n,k,k') = l/(27r) 



(11) 
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with the assumption that both k and k' lie within the first 
Brillouin zone so that their difference is not a reciprocal 
lattice vector jl8j |. 

The case n ^ n' can be carried out by considering the 
eigenvalue equation for the w„(k,x) functions 

^o(k)|« Tlik )=^(k)|« n|k > (12) 
where we have defined Ho(k) as (p = —ifi'V) 

H (k) = ^-(p + hk) 2 + V L (x) (13) 
2m 

Then, by left multiplying Eq. (fT2")l by (u , k | and using 
the equivalence 

ffo(k) = H Qs!) + ^-{k 2 - k 12 ) + -p-(k-k') (14) 
we get 

with the momentum matrix elements P nn >(k, k ) defined 
by 

p nn ,(k,k') = ^- / <(k,x)(-zWK/(k',x) 

" J cell 

(16) 

and 

A£w(k,k') = E n (k) - E n ,(k') - |- {k 2 -k' 2 ) (17) 

2to v ' 

Finally, the equation (JT]) for the expansion coefficients 
can be rewritten as 

ihd t ip n (k) = E n (k)cp n (k) + [ U(k-k')ip n (k') (18) 

JW 

where it is easy to identify the single band dynamics 
(first line) and the interband coupling (second line). This 
equation is so far very general and relies on the only as- 
sumption that the external potential U as no appreciable 
variation on the scale of a single lattice cell (see Eq. ©). 

Let us now transform back the above equation in co- 
ordinate space; this can be achieved by projection on the 
Wannier basis 

*( x ) = EEx«(^vr(x-R«) (19) 

where the Wannier basis functions satisfy the orthogo- 
nality relation 

f 0jf*(x-R i )^(x-R i ) = 5 nn ,5. l3 (20) 

J X 



and can be expressed in terms of Bloch functions as 



^(x-R i ) = y^3^Vn(k,x-R i ). (21) 

The use of the Wannier basis has two advantages: (i) 
the amplitudes Xn(Ri): that play the role of envelope 
functions on the new basis (see Eq. (fT§|) ). can be ob- 
tained from the Bloch coefficients in Eq. (p~8|) by a simple 
Fourier transform 



X " (Rl) = V(2^/^" (k)eik ' Rl; (22) 

(ii) they can be interpreted as the actual wave function of 
an electron in the nth band if one is interested in "macro- 
scopic" properties of the system on a scale much larger 
that the lattice spacing (that is equivalent to average on 
a scale of the order of the lattice cell). For example, by 
using the completeness of the Wannier basis in Eq. (|20p . 
the density and current distributions can be expressed as 

n 

J, = (J(x)) ce «_i ~ Ai m V[ X :(Ri)Vx„(R i )](24) 
im £ — ' 

n 

Since the functions Xn(Rj) are in principle defined only at 
the lattice sites, it is convenient to follow the approach of 
(25| and perform the limit to the continuum by extending 
the dependence of the Xn(Ri) to the whole space (Rj — ► 
x). This yields the following expressions for the cell- 
averaged charge and current densities 

p(x) ~ ]T|x„(x)| 2 (25) 
J(x) ~ Ai m ^[ x;(x)VXn (x)] (26) 

n 

Then, by using standard properties of the Fourier 
transform, Eq. (fT8"j) can be formally written in coordinate 
space as 

ihdtXni*) = E n (-ihV)xn(x) + Z7( x )Xn(x) 

* nr r . kx r p nn ,(k,kQ 

m^J (2„)»J k J k , AE nn ,(k,k') 

.(k-k')L>(k-k')^(k') (27) 

This equation is equivalent to the generalized form of 
the Wannier equations of Ref. [25| , with the advantage 
of having the interband term written in a more trans- 
parent form in terms of its Fourier components. This 
expression allows for a simple manipulation of the above 
equation, that for practical use, has to be further sim- 
plified. The simplest approach is to adopt the following 
standard approximations [26j , assuming that 
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i) the energy spectrum is of simple form with min- 
ima/maxima of each band at some point k = k in 
the first Brillouin zone; 

ii) the (f n (k) functions are localized on a small region 
of k space around k — fco during the whole evolu- 
tion of the system. 

For convenience in the notations, and without loss of gen- 
erality, in the rest of the paper we will set fco = 0. Then, 
let us consider the term P nn ' (k, k )/ AE nn ' (k, k ) that 
characterizes the interband coupling; to first order in k, 
k' we can write 



,(k,k') 



AE nn ,(k, k' 
1 



AE„ 



) AE nn > 
k-V fc P n „/ 



(28) 



k'-VvP, 



0(k 2 ) 



where to simplify the notation we have defined AE nn / = 
A£„„,(0,0) = E n (0) - E n ,(0) = E n -E n ,, P nn , = 
P nn /(0,0), and we have used the fact that the energies 
do not contain first order terms in k (see (ii)). 

The first derivatives in Eq. (|28|) can be evaluated by 
using the relation 

i 



that can be obtained by differentiation Eq. (fl2|) with 
respect to k, a nd p rojecting the term Vfeti n (k,x) on the 
u n (k,x) basis [23]. Then a straightforward calculation 
yields 

VfcPw = A J2 P r" P r"' - Mnn< (30) 



Vfc/P r 



171 IT! E n~ E n" 

n" 

Tl x "s -Pnn" Pn"n' 



E 



m j-f E n >-E n „ 



M* >n (31) 



where P*, n (k,k') = P nn /(k,k'). The last term of Eq. 
(|27p then becomes 



Pnn' (k, k 



k' AE nn ,{k,\sl 

nn 



{-•(k-k')L>(k-k')^„'(k') 



(32) 



I E ir- I (k-k')c/(k- k ')^(k') 



n'^n 



nn' Jk 



+ ^ E ^ / (k-k')f/(k-k')kVn<(k') 

n'^n 

+k- E ^ / (k-k'Mk-k' W (k') + o(k 2 ) 
This expression allows us to write Eq. (|2T|) as 



ift&X„(x) = £ n (-iW)x„(x) + (7(x) X „(x) - iV£/(x) - V P "" 

77? ^ 



/a " AE nn 

n'^n 



■Xn'(x) 



(33) 



E if^ V ^(x)- A ^ ^l[V 2 t/(x) X „,(x) + VC/(x)Vx„(x)] 

/ft- — 1 JJj<ri ti 1 lib A Pjin t-t 1 



m *—f AE nn 



m '—f AEnn 1 

n =£n 



The above equation represents the main result of this 
paper. It describes the evolution of the Wannier envelope 
functions by fully including the effects of the periodic 
potential and accounting for the interband coupling due 
to the perturbation potential U up to second order in k. 
Eq. (|33|) can be further simplified by means of the usual 
effective mass approximation that amounts to retaining 
only up to quadratic terms in k in the kinetic operator. In 
general this corresponds to replace the bare mass by a 3 x 
3 mass tensor m*j [23| ; in the special case of an isotropic 
periodic potential or for a one-dimensional system as we 
will consider later, one simply gets 



E n (k)=E n + ^+0(k 3 ). (34) 



III. COMPARISON WITH OTHER k P 
MODELS 

Let us now discuss the connection of the model in Eq. 
(|33| with other two "k-P" models, namely the Kane and 
Luttinger-Kohn (LK) models [HI, H2] , which are widely 
used both to estimate the band diagram in semiconductor 
and the transmission coefficients of interband devices [l?], 

m. 

The Kane model is based on the following choice of the 
basis elements 

(x|n,k)jf a = e lk ' x u n (0,x) , (35) 

that form a complete orthonormal set, and can be used 
to expand the electron wave function f in a similar way 
of what shown in the previous section. On this basis the 
Schrodinger equation takes the form 

ihdtfn (k) = E / H ™> (k, k') <p n . (k') (36) 
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where the hamiltonian matrix elements are 
Hf„ Q ,(k,k') = (n,k\H Q + U\n',k') Ka 



(37) 





h 2 k 2 \ 






2m / 



m 



5(k-k') 



-17(k-k')<J n 



By means of an inverse Fourier transform (see Eq. 
([22]) ) it is straightforward to recover the equation for 
Kane envelope functions x^" a (x) 



2m 



V^ + £„ + C/(x) x r(x) 



-iA^ Pw .Vx^°(x). 



(38) 



This equation shows that in the Kane representation en- 
velope functions related to different "band" indices are 
coupled even if the external field is vanishing. This is due 
to the fact that the unperturbed hamiltonian H is not 
diagonal on the Kane basis (see Eq. (37])), and therefore 
the n here does not correspond to the usual band index 
of the Bloch picture. In other words this means that the 
envelope functions x^" a ( x ) do not have the direct physi- 
cal meaning of wavcfunctions of an electron in a definite 
energy band. As a consequence, one should be careful 



in estimating truncation errors when the full problem is 
reduced to a finite set of envelope functions. 

To overcome the previous difficulty Luttinger and 
Kohn proposed a different choice of the basis functions 
fl9| . The idea is to use a quasi-unitary transformation 
to diagonalize the Kane hamiltonian in the momentum 
space up to first order in k. In this way, it is possible to 
get a natural extension of the effective mass single band 
model in the multiband framework. The new hamilto- 
nian reads 



n LK = Q- 1 n Ka e 

where is defined as follows 



(39) 



(n,k|6|n',k')Ka 



m AE nr 



(40) 

providing a unitary transformation to first order in k. 
Accordingly, the elements of the LK basis are defined by 
\n, \s)lk = 0| n ,k)xa, and correspond to an expansion 
of the u ra (k, x) functions to first order in k 



(x\n,k) 



LK 



ikx 



u„(0,x) + k 



<9u„(0,x) 



<9k 



(41) 



In the coordinate space the LK model reads 



iw tX i K {x) = 



E n — 



2ml 



U(x) 



LK 



m *—f AE nn > 
n / n 



(42) 



— V (PwV) (P n »n'V) 

711 n ^ ' 



TO ,,, 



AE„ 



X^(x) 



The first line here corresponds to the first line of Eq. |33]) 
with the effective mass approximation (]34]k the second 
line instead represent a spurious coupling between differ- 
ent bands that corresponds to the choice of the truncated 
basis in Eq. (j4i"]) . and is usually neglected [HI]. In our 
approach this term would come from the expansion of 
the off-diagonal kernel 

(n,k|if |n',k')= / e- ik - x <(k,x)if e ik x W „/(k',x) 

J X 

but is cancelled exactly by a term coming from the ex- 
pansion of the u„(k, x) functions to second order in k. 
As a matter of fact, this contribution is absent in our 
approach since only the off-diagonal terms that depend 
on the external potential U have been approximated (up 
to 0(k 2 ) in Eq. (]33]l). 

We also remark that in the LK approach one usually 



neglect also the interband coupling proportional to the 
applied field VJ7, and this prevents any description of 
interband tunneling effects. 



IV. AN EXAMPLE 

As an application of the model discussed in Sec. HTlwe 
consider a one-dimensional semiconductor device consist- 
ing of a multilayer heterostructure where only two bands 
play a relevant role, namely the "conduction" and "va- 
lence" bands. As a further approximation we keep the 
interband terms only to first order in k, neglecting the 
terms proportional to the matrix M nn /, and adopt the 
effective mass approximation (1341) . Thus the system in 
Eq. (]33j) can be reduced to the following set of coupled 
equations 
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ihdtXc(x) = E cXc {x) - ^—^ 2 Xrix) + U(x)xc{x) - iVU(x)^rXv(x) 



2m 



HP 



(44) 



ihdtXv(x) = E v Xv(x) + -, — -jV Xv(x) + U(x)xv(x) - iVU(x) — —Xc{x) 



2 m 



mE n 



that depend on four phenomenological parameters: the 
interband momentum matrix P = P c v = P* c (see (Tol . 
[24l | for a numerical estimate) , the energy gap E g = E c — 
E v , and the effective masses m* v for the conduction and 
valence bands respectively (m* = — |m*|). 

The total potential can be written as U(x) = Uh(x) + 
U e (x), where U e (x) is the electrostatic potential gener- 
ated by the charge distribution in the device, and Uh(x) 
accounts for the spatial dependence of the band edges in 
the heterostructure. Indeed, in real heterostructures the 
conduction and valence band edges depend on x and to 
account for this we adopt the point of view of considering 
the Bloch spectrum as constant among the layers, treat- 
ing the actual spatial dependence as an external poten- 
tial applied to the heterostructure bulk. We remark also 
that here we are considering only coherent transport, ne- 
glecting any dissipative phenomena like electron-phonon 

I 



scattering that are not expected to affect significantly the 
tunneling process. 

The electrostatic potential can be calculated self- 
consistently by using the Poisson equation 

eV 2 U e (x) = qp{x) = q 2 [C{x) - n{x)] (45) 

where the total charge distribution p(x) is the sum of the 
charge concentration qC(x) of the doping ions plus the 
charge distribution —qn(x) of free electrons. 

Let us now consider an heterostructure device in con- 
tact with a source and a drain reservoir at a temperature 
T. In the presence of an incident electron beam with mo- 
mentum q and energy E{q) injected from the reservoirs 
into the device, the steady state of the system is obtained 
from the solution of the stationary equations 



E(q)xl(x) = E cX l{x) - A-V 2 X ?(.t) + U{x)y*{x) - iVU(x)^ X q v (x) 

(46) 

E(q)xUx) = E vX l{x) + ^—W\l{x) + U(x)xl(x) - iVU(x)-^rxi(x) 
I 2 m* mE a 



combined with the Poisson equation (|45|) . The equations 
are solved by approximating the spatial derivative by 
a Runge-Kutta method and using a Gummel predictor 
scheme to reach convergence [30} . 

The electronic density n — n c + n v is constructed in 
terms of the pure state solutions x% an d xl 01 t ne above 
equations, weighted by the momentum distribution of the 
incident beams 



n(x)= / dqf (q) [\ X 9 c (x)\ 2 + \xl(x)\ 



(47) 



where fo (<?) is the Fermi-Dirac distribution integrated on 
the transverse coordinates [3lj ] - 

Similarly the electronic current J = J c + J v is calcu- 
lated as 



dq fo(q)Imlxl(x)S7xUx)} . (48) 



To model the charge injected in the device from the 
source and drain reservoirs we use transparent boundary 
conditions. For example, at x = 0, in case of electron 
beam incident in the conduction band with positive mo- 
mentum q, we have 



d 

dx 
were 



x=0 



iki 













^[E(q)-E v ] ; 



(49) 

(50) 
(51) 



the other cases are treated in similar way. 

Let us now discuss what happens in case of that is 
a one-dimensional RITD, where charge transport across 
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FIG. 1: (Color online) Simulated heterostructure profile and 
doped regions of the RITD. The widths of the layers are cho- 
sen: 5.00 nm for the quantum well (Q.W.), 3.00 nm for the 
barriers (Bar.) and 4.50 nm for the spacer layers (Spac). 
Resonant tunneling takes place when the energy of electrons 
in the conduction band is resonant with that of the bounded 
hole state in the central quantum well (shaded area). 



the device takes place thanks to resonant tunneling be- 
tween electron states in the conduction band and hole 
states in the valence band. As a specific case here we 
consider the simple test device depicted in Fig. [TJ con- 
sisting of a 5.00-nm-wide quantum well (for hole states), 
bounded by two identical 3.00-nm-wide barriers that are 
interfaced to the bulk by two 4.50-nm-wide spacer lay- 
ers. In Fig. [1] we also show the resulting heterostructure 
potential in case of a GaSb lattice bulk with a doping 
concentration of 10 18 cm~ 3 32]. For consistency with 
the formalism presented in the previous sections, we are 
considering a constant band gaps heterostructure (we re- 
mind that in our model all the interband effects are taken 
in account by an external potential added to the bulk 
Hamiltonian) . We note that, as far as resonant tunnel- 
ing is concerned, the actual shape of the conduction band 
far above the resonance level cannot produce strong ef- 
fects on the tunneling process. Therefore we expect that 
the I-V behaviour does not change qualitatively in case 
of more complex conduction bands profiles. 

In Fig. [2] we show the calculated equilibrium self- 
consistent potential U and the density of electrons cor- 
responding to the unbiased case: note that, in this case, 
the profile of the heterostructure potential Uh of Fig. [T] is 
practically unchanged by the addition of the electrostatic 
potential U e . 

The steady TV characteristic of the device at a tem- 
perature of 300 °K is shown in Fig. [3l where the current 
I flowing through the device is plotted as a function of 
the bias voltage V b applied to the drain contact. This 
picture shows that the model is capable to reproduce the 
expected negative differential resistance (NDR) in a cer- 
tain range of values of the applied potential (here for 




x (nm) 

FIG. 2: (Color online) Self-consistent potential profile and 
density of electrons corresponding to unbiased case. 




FIG. 3: (Color online) I — V characteristic of the simulated 
diode. Notice the negative differential resistance for Vb > 
V = 0.05 V. 



V b > V = 0.225 V). 

In Fig. 0] we show the calculated self-consistent po- 
tential profile U and the density of electrons correspond- 
ing to peak and valley currents, for Vb = 0.225 V and 
Vb = 0.27 V respectively. We note that as expected the 
density of electrons in the central region of the potential 
profile (see Fig. 0]) is an increasing function of the ap- 
plied bias below resonance (V& < Vq), and then sharply 
decreases in the NDR region. 

These results demonstrate that the present model is 
promising for reproducing the TV curves of physical de- 
vices as that considered e.g. in [H, [Hj], although a pre- 
cise comparison requires the extension to structures with 
varying band gap profiles (or, more generally, with x- 
dependent Luttinger-Kohn parameters [2(|) and will be 
considered in a future publication. 
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FIG. 4: (Color online) Self-consistent potential profile and 
density of electrons corresponding to the peak (top) and to 
the valley (bottom) currents. 



the usual Bloch theory by means of a /c-expansion, and is 
formulated in terms of cell-averaged envelope functions 
obtained by projection in the Wannier representation. 
The model is suited to describe in a clear fashion tun- 
neling effects between different bands in presence of an 
applied potential. Its advantages with respect to other 
widely used approaches has been discussed. 

As an application we have considered the case of a 
RITD, a heterostructure device where the electronic cur- 
rent flows between a "conduction" and a "valence" band, 
interfaced by potential barriers. In this case the model 
is reduced to a system of two Schrodinger equations for 
the electron envelope function coupled with the Poisson 
equation for the field generated by the electronic distribu- 
tion itself. It nicely reproduces the expected behaviour of 
the current as a function of the applied voltage, exhibit- 
ing a negative differential resistance in a certain range of 
values of the applied bias. 

The present approach is therefore particularly promis- 
ing for reproducing the behavior of physical devices char- 
acterized by resonant tunneling, and we are currently 
working at the extension of the model which takes in 
account the specific symmetry properties of the crystal 
lattice and to treat degenerate and varying band gap pro- 
files. Moreover, with the inclusion of spin-orbit coupling, 
our model may also be relevant for the field of spintronics, 
where k-P methods are becoming increasingly important 
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